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PREFACE 


This work was undertaken because an analytical method used to predict 
lunar and solar perturbations for earth satellites failed to give suitable values 
for the perturbations in the orbital elements of the satellite Relay II. A study of 
this problem has yielded a feasible method by which one may accurately compute 
perturbations in the angular elements for Relay II. This study partially satisfies 
the requirement toward a Masters Degree at Georgetown University. 

I wish to express my especial gratitude to Mr. David Fisher of Goddard 
Space Flight Center for the guidance he has given me throughout this work. In 
addition, I am indebted to Father Francis Heyden, S. J. of Georgetown University 
for the unfailing interest he has shown in this work and in all my studies at the 
University. I also wish to thank Dr. Thomas Margrave of Georgetown University 
and Mr. Theodore Felsentreger of Goddard Space Flight Center for their advice 
and encouragement. 
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INTRODUCTION 

Six orbital elements, coordinates, at any instant completely specify an orbit. 
These six elements or coordinates are usually the three Cartesian coordinates 
of pcsition and the three of velocity or the six Keplerian elements. Since earth 
satellites can in the first approximation be assumed to move in elliptical orbits, 
the Keplerian elements form a suitable set of coordinates with which to describe 
a satellite's motion. In a simplified case, the two-body problem, there is a com- 
plete analytical solution. In other words , there is a solution given by a simple 
algebraic formula which yields the value of the six elements at any desired time. 
The actual orbit problem for a satellite, however, is much more complicated. 

The solution tc the problem of describing an actual satellite's orbit is approached 
either numerically or by use of an infinite series technique. The use of such 
methods still leaves the Keplerian elements as suitable quantities with which to 
solve problems but their use can occasionally cause difficulties. 

One such difficulty arises when one attempts to describe the motion of the 
artificial satellite Relay II. The period of an argument of the motion, the longi- 
tude of perigee, is about 570 years. This condition is referred to as resonance. 
Resonance complicates the computations of the Keplerian angular elements u> and 
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ft for Relay II. The series expressions for u (or g which is the argument of 
perigee) and ft (or h which is the longitude of the ascending node) have some luni- 
solar terms which have the derivative of the longitude of perigee in the denomi- 
rator. Since the longitude of perigee for Relay II has a value of almost zero, 
such luni-solar terms produce a small divisor problem. These near singular 
terms create sizable errors in the computed values for the orbital elements. 

This work shows the feasibility of an approach to treat the angular orbital 
elements u> and ft for Relay II. The approach numerically integrates Lagrange's 
planetary equations and thereby bypasses forming terms which have the deriva- 
tive of the longitude of perigee in the denominator. By using this method, one 
avoids all near singular terms which cause large errors. The approach has 
been checked by two computer programs and has been compared to the work 
of Theodord Felsentrager. A determination of the J 3 zonal harmonic of the 
earth's gravitational potential was made as a check to the work. 


PART I 


DERIVATION OF THE EQUATIONS DESCRIBING THE MOTION OF 

THE SATELLITE RELAY U 


CHAPTER I 

FORMULATION OF THE TWO-BODY PROBLEM 

Due to the rapid development of satellite theory, the problem of describing 
the motion of a satellite which has a near-resonant orbit condition seems rather 
removed from basic celestial mechanics. The equations which are finally arrived 
at to describe the behavior of the angular orbital elements for Relay II are 
actually simple developments from such fundamental laws as those of Newton. 

A detailed development starting from these basic laws will give physical insight 
to the problem and acquaint the reader with the notation. 

Newton's first law states that if the resultant force acting on a particle is 
zero, the particle will move with constant velocity. Where this law holds, we 
have an inertial frame of reference. The first law is described by the equation 

T = 7 0 + f t (1) 

where r is constant. The second law states that if the resultant force acting on a 
particle is not zero, the particle will suffer an acceleration proportional to the 
magnitude of the resultant force and in the direction of the resultant force. In 
vector notation 
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• • 


(2) 
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= d(m r) - 
dt 


The third law states that when two bodies act on each other, the forces they exert 
will be of equal magnitude but opposite in direction. This is described by the 
equation 


F 2 i - F 12 . (3) 

In addition to these three laws, there is Newton’s law of gravitation which states 
that any two particles of mass m 1 and m 2 separated by a distance ? 12 mutually 
attract each other with a force directly proportional to the product of their masses 
and inversely proportional to the square of the distance. The force exerted on 
particle one by particle two can be written as 


k 2 


F 12 = 


m i m 2 r 1 2 

F 13 
1 12 1 


(4) 


By using the previous four equations , one can immediately investigate the 
simplest orbit problems. One such problem can be described as occurring when 
two bodies, each of whose mass is spherically distributed, orbit about a mutual 
center of gravity. Let these masses be denoted and m 2 . The external gravity 
field for a body with a spherical mass distribution is the same as for a body 
having »\11 its mass concentrated at its geometric point center. 

The center of mass of the two-body system is defined as 


2 



i = 1 


m. 

1 
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If the origin of the coordinate system is taken to be at the system’s center of 
mass, then the motion of m 2 relative to m 1 can be obtained by setting equations 
(2) and (4) equal, i.e., 


and 


See Figure 1. Let 


~ k 2 m l m 2 ( p 2 - p x ) 
m i P\ - 

l p 2 -pJ 3 


- k 2 mj m 2 (p 2 - pj 



r 


1 2 


P 2 " Pj 


Then 


- k 2 (nij + m 2 ) r 


12 


12 


r 12 l 3 


(6) 


( 7 ) 


(8) 


If mj and m 2 can both be considered point masses, then equation (8) has a 
solution 


a (l - e 2 ) 
r 12 1 + e cos f 


O) 


This equation describes a conic section; see figure 2.* Only the elliptic orbit 
case will be considered here (0 < e < 1). In this case, three orbital elements 
can be immediately defined. The size of the ellipse is given by the value of its 


forest Moulton, An Introduction to Celestial Mechanics (New York; Mac- 
millan, Co. , 1935), p. 321. 




APIFOCAL DISTANCE = a (1 


- 

+ e ) X 


PERIFOCAL DISTANCE ^ a (1 - e) 


Figure 2. — Elliptical orbit of a satellite with respect to a primary body 
located at focus f,. 
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semimajor axis, a . The semimajor axis is proportional to the total energy 
in the orbit. The shape of the ellipse is described by the value of the eccentricity, 
e. From figure 2 

f* ~ c/a = j/a 2 - b 2 /a . 

The quantity f is termed the true anomaly. It is the angle formed at the focus 
f j between the radius vector of the satellite at any time t and that to the peri- 
focal point at time T, which is called the time of perifocal passage. These 
definitions can be used to calculate the perigee and apogee heights of the satel- 
lite Relay II (a - 1.7449 r e and e = .239). The 

Apogee height = (1.7449 r e (1 + .239) - r e ) (3963 mi./r e ) 

= 4728 miles 

while the 

Perigee height = (1.7449 r e (1 - .239) - r e ) (3963 mi./ r e ) 

= 1299 miles. 

For the two- body problem, a complete analytical solution can give the position 
of the satellite at any specific time in the past or present. It is not necessary 
to do a step by step time integration to obtain the coordinates at a desired time 
(which is one method of solving problems involving more than two bodies). For 
the two-body problem, any of several sets of six independent coordinates will 
suitably describe the orbit. For example, the six coordinates x, y, z, x, y, z 
at any one time would specify an orbit. The two-body problem however has a 
preferred set of elements, the Keplerian elements. They are convenient because 
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the time derivatives of the Keplerian elements are then zero. In other words, 
when there are no perturbations, the Keplerian elements are constant. They 
thus describe two-body motion for all time.* 

The Keplerian elements a and e which specify the size and shape of the 
ellipse have been mentioned. To define the others, refer to the figure 3 which 
depicts a celestial sphere. The projection of the earth's equator on the celestial 
sphere forms one fundamental great circle, the celestial equator. The projection 
of the sun's apparent path during the course of the year defines the second funda- 
mental great circle, the ecliptic. The point where the sun moving from south to 
north along the ecliptic crosses the equator is called the vernal equinox, y . 

A suitable coordinate system can be set up by pointing x toward y and z toward 
the north celestial pole. The satellite orbit olane as seen from the center of the 
earth projected against the sky is also plotted. The spherical angle between it 
and the celestial equator defines another Keplerian element, the inclination, i. 

The angular distance measured westward along the equator between y and the 
point where the satellite crosses the equator traveling in a northerly direction 
defines the element 0, the longitude of the ascending node. The angle between the 
ascending node and the perigee point on the orbit is referred to as a;, the argument 
of perigee. Any of several elements can be taken to be the sixth Kepler element. 

If the mean motion of a satellite , n , is given then the sixth orbital element 
generally used is the mean anomaly, M. The mean motion, 

*Wolf Research and Development Corporation, "Orbital Dynamics," Satellite 
Geodesy, Theory and Applications (College Park, Maryland: Wolf Research and 
Development Corporation, 1968) VI, p. 1. 



Figure 3. -Projection of the satellite's orbit upon the celestial sphere. 



11 


n 2~ P (10) 

where P is the orbital period. The quantity M is defined as the product of the 
mean motion and the time elapsed since the perigee passage, ( t - T). A hybrid 
element F , the longitude of perigee, is defined as the sum of 0 and oj. 

Any position vector F and velocity vector t at time t has a corresponding 
unique set of Kepler elements a, e , i , M, and H. For any instant, there is a 
coordinate transform relating the two sets of coordinates. The problem can be 
viewed in either of the two sets of coordinates. 


CHAPTER II 

MODIFICATIONS TO THE TWO- BODY PROBLEM IN ORDER TO 
DESCRIBE THE MOTION OF A REAL SATELLITE 

The two-body problem does not describe the actual motion of a satellite 
like Relay II very well. This is because the earth does not have a spherical 
mass distribution. Consequently the gravity field of the earth cannot be con- 
strued as acting from a point source located at the earth’s center. Also the 
sun and the moon significantly effect the orbit of Relay II. For a close earth 
satellite, the simple two- body problem still remains a good first approximation 
to the actual physical solution. Since these two solutions are close for Relay II , 
the actual solution can be considered to be made up of a two-body motion and a 
deviation from the two-body motion, a perturbation. It is inconvenient to describe 
this motion analytically in terms of the x, y , z , x, v , z coordinates; they change 
too rapidly. The Keplerian elements a, e, i, M, fi on the other hand change 
only relatively slowly. 1 At any instant, there exists a set of Keplerian elements 
that would be followed if there were no forces other than the two-body force. 


^bid., ’’Orbital Dynamics,” VI, p. 1. 
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These are called osculating elements. When disturbing forces do act, the co- 
ordinates in the actual orbit are changed from those which the body would have if 
there were no disturbing force. This change in position is referred to as a per- 
turbation in the coordinate. A perturbation (as S g ) must be added to the two-body 
value of the coordinate (as g ) to obtain the actual value g + Sg. 

The analytical expression for the motion of Relay II can be obtained by 
assuming that the perturbations to the two-body motion come from three sources. 
The equations describing the perturbations due to two of these sources, the sun 
and the moon, can be derived by considering a simple form of the three-body 
problem. The perturbations due to the third source, the nonspherical earth, must 
be considered separately. In the simple version of the three body problem, we 
assume that there are two sizable masses, and m 3 ; each has a spherical mass 
distribution. Let a third body m 2 , of negligible mass, orbit around body nij. Let 
a remote body m 3 perturb the two* body motion of m 2 about m 1# To obtain the 
equation of motion for m 2 , the forces given by Newton’s second law and the law 
of gravitation are equated. The total force on body one is 


„ k 2 

m l m 2 _ k 2 

m l m 3 

F, = m, r, = p 12 +— — p 


•l *i 


Pl2 ]3 


0 13 

^13 


13 • 


( 11 ) 


The total force on body two is 


F 2 - m 2 r 2 - - 


k 2 nij m 2 _ 

l^ 12 l 3 


k 2 m 2 m 3 _ 


12 


'23 


P 


23 


3 


( 12 ) 
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By transferring the origin of coordinates to nij , one obtains the equation of motion 
of m 2 relative to m x , i.e., 


P 12 



f* 


12 


k 2 


m. 


Tj = - k 2 (mj + m 2 ) ' -f rw iii^ 

lp, 2 l 3 




(13) 


where 


k 2 (mj + m 2 ) = /lx . 


Equation (13), describes the motion of m 2 as being composed of a two-body motion 
plus a perturbation. 1 Using the notation of figure 5, let 

p i2 = r , 

P 13 = r 

and 

p 23 = A. 

In terms of these variables 



(14) 


Let R denote the disturbing function and U denote the potential; then the above 
equation can be written as 


^loid., ’’Coordinate Reference Systems Used in Astrodynamics,” II, p. 7. 
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Figure 5. — The three-body problem consisting of the earth, the satellite 
and the moon. 
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k 2 (mj *m 2 )r 


lr 1 3 


Integrating equation (15) yields the forcing function (or the potential) 


k 2 (m 1 + m 2 ) 

U r r-r + R • 


The disturbing function 


k 2 m„ I r=r - 


1 r ' • r 


3 VIA! T' 3 


where 


r cos i r 


The term (k 2 m 3 /A) is called the direct perturbation and the term [k 2 m 3 (|7|/ 

| 7' j 2 ) cos 77'] is called the indirect perturbation. 

Henceforth the vector notation will be dropped so as not to complicate notation. 
To obtain an expression describing the effect of luni-solar forces on a satellite, 
one must visualize the following situation. A perturbing body, mass m 3 , acts 
upon an earth satellite, mass m 2 . The perturbing body m 3 can be either the sun, 
mass m r , or the moon, mass m c . The direct perturbation on the satellite due to 


m 3 1S 




X S. W. McCuskey, Introduction to Celestial Mechanics (Reading, Massa- 


chusetts: Addison-Wesley Co., 1963), p. 101. 




IS 


By the law of cosines, 


:. 2 = r 2 


2r r ' 


cos 



then 


|_1 “ 2 ~v T cos * (tO J 

This equation can be expanded into an infinite series 

c c 

3 = 1 - 2^ (r) Pi (cos rr' ) . (18) 

i = 1 

The term P 4 is the Legendre polynomial of the i th order. The disturbing 
function then becomes 


m. 


■L 

i = 1 


(f) 1 p.( 


cos rr 




cos rr 


(19) 


There is no simple analytical solution describing the motion of a body in a 
three or more body system. Numerical and infinite series approaches must then 
be used. The use of equation (19) for R is an infinite series approach. Kozai 


*D. Brouwer and G. Clemence, Methods in Celestial Mechanics (New York: 
Academic Press, 1959), p. 119. 

2\\ T olf Research and Development Corporation, ’’Orbital Dynamics,” X, p. 6. 
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used this disturbing function to develop a formula describing luni-soiar effects 
to the order i = 2. 1 Murphy and Felsentreger’s paper extend this to i = 3 but 
they consider only long-period terms (terms independent of the mean anomaly). 2 
Among these terms are the two long- period luni- solar resonant terms for Relay 
II. These two terms are given in Felsentreger’s paper as F R . 3 They are eferred 
to in this work as R R(SM) . The resonant part of the luni-solar disturbing func- 
tion is 

R r(SM) = + 64’ a2e2 ( 1 " 2 cos i T cos2 OjVi 2 m s sin 2 i© cos 2(^7 ~ g - h ) 


+ 



1 

2 sin 


cos 2( g +h) 


- sin i" cos i sin i,.(l _ cos i^cos (k" + 2g + 2h) 
* sin i" cos i" sin i c (l+cos i„)cos (H - 2g _ 2h) 


+ ^ sin 2 i " (l cos i z ) 2 cos 2(ft" + g + h) 


^ sin 2 i (l + cos i^.) 2 cos 2(Cl" - 


g- h )Jj- 

( 20 ) 


ly. Kozai, ”On the Effects of the Sun and the Moon Upon the Motion of a Close 
Earth Satellite,” SAP Special Report, Number 22 (March, 1959). 

2 J. P. Murphy and T. L. Felsentreger, ”An Analysis of Lunar and Solar 
Effects on the Motion of Close Earth Satellites,” NASA TN-D-3559 , (August, 1966). 

3t. L. Felsentreger, ’’Long Period Lunar and Solar Effects on the Motion 
of Relay n ” GSFC X-547-66-102, (March, 1966), 3. 
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To obtain simplified solutions in the previous cases, it was assumed that 
each mass was spherically symmetric. In such a case, the body would attract 
an external body as if all its mass were concentrated at its center. To repre- 
sent reality more closely one must complicate such a simple mass model. An 
accurate mass model must take into account any mass distributions which cause 
the gravitation field of the earth to significantly deviate from a spherical po- 
tential. One model for the earth’s mass separates the irregular mass distribu- 
tion which deviates from a spherical distribution into mass components that are 
symmetric to the axis of rotation. The major component in this model is the 
equatorial bulge of the earth (the earth’s equatorial diameter exceeds its polar 
dimaeter by some 27 miles). The earth also has a mass distribution component 
deviating from sphericity which is referred to as pear-shaped. This small 
effect only amounts to a 50-foot departure at maximum from a spherical earth. 
See figure 6. There are many other significant mass distribution components. 
Such axially symmetric mass distribution components, which account for the 
earth’s irregular mass distribution, produce the various zonal harmonics in the 
earth's gravitational potential. * 

The equation of motion for a satellite moving in a nonspherical gravitational 
field of potential U is 

r = V U . (21) 

*C. A. Wagner, ’’The Drift of a 24-Hour Equatorial Satellite Due to an Earth 
Gravity Field Through 4th Order.” NASA TN-D-2103 , (February, 1964). 


21 


J 2 ZONAL HARMONIC 


J 3 ZONAL HARMONIC 


J 4 ZONAL HARMONIC 


NP 

t 

4.3 MILES 


2.1 MILES 


NP 

t 



t 




Figure 6.— The primary mass distribution components symmetric to the 
axis of rotation which depart from a spherical earth. 
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The nonspherical potential can be written as 

J = - + R. (22) 

r 

The term (^/r) is the potential due to a spherical earth and R is the contribu- 
tion to the potential arising from mass distributions that deviate from a spherical 
distribution. 

To obtain an expression for the potential of an irregularly shaped body, one 
resorts to Laplace's equation. The potential must satisfy Laplace's equation 
outside the irregular body, the earth, i.e., 

V 2 U = 0. (23) 

In spherical coordinates, equation (23) becomes 


1 d ( , dlA 

r ~7Tr W Jrj 


1 a 2 U 
r 2 sin 2 - d : 2 


1 

r 2 sin 



^sin 



0. 


If one considers zonal harmonics only, then the mass distributions are symmetric 
to the axis of rotation and U is independent of i-. The general form of the po- 
tential U is then obtained by solving Laplace's equation assuming a solution of 
the form 

U (r, 6) = R (r) 0 (5) . 


A convenient form of the final solution is 


U = 


m k 

earth 


2 r r 2 r 3 

1 - - 1 J 2 P 2 (cos 6) - — J 3 P 3 (cos 6) ~ • • 
r 2 r 3 


(24) 
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where P 2 is the Legendre polynomial of the second order. 1 The quantity 


P 2 (cos 6) = — (3 cos 2 Q - 1) . 


The quantities J 2 , J 3 and J 4 are the spherical harmonic coefficients of the 
earth’s gravitational field; their magnitudes indicate the strength of the com- 
ponent of deviation from a spherical gravitational field. The first term in the 
brackets is due to a spherical earth; the second is due to the equatorial bulge and 
the third term, the J 3 term, is due to the pear-shaped earth. 

Using this above expression of the potential, Brouwer derived expressions 
to suitably describe the motions of the Keplerian elements for earth satellites. 


= 90 - 8 


m e a r t h k2 


Then Brouwer’s potential function takes the form 


U = — 


Jw 


u V ™ 1 J k 1 « 

- '-L — 


P. (sin 8) 


By spherical trigonometry, 


sin 8 = sin i sin (g + f); 


J J. M. Danby, Fundamentals of Celestial Mechanics (New York: Macmillan 
Co., 1962), pp. 108-111. 

2 Brouwer and Clemence, Methods in Celestial Mechanics, p. 563. 
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tl*en 


fJ- 2 r e 17 1 3 2 A a 3 

u = t + ^3 [y 3 + 3 cos2 v 7^ 


cos 2 ■— cos (2g + 2f)j 


+ • ' 


(26) 


In addition to this formula, Brouwer needed Hamilton's equations 


dg *i _ dF 

dt “ <9p i 

and 

d Pj _d¥_ 

dt dq { 


(27) 


(28) 


These two equations determine the changes dq £ and dp. in the coordinates q. 
(x, y , z ) and in the momenta p. (mx ,my , mz) that occur during an infinitesimal 
period of time dt. In these equations F is the constructed Hamiltonian of the 
system. It is equal to the negative of the total energy of the system, i.e., 

F( Pl . q,) = U - T ke = £ + R • (29) 


*D. Brouwer, "Solution of the Problem of Artificial Satellite Theory Without 
Drag," Astronomical Journal, LXIV, p. 378. 
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In the orbit problem, the coordir.ates and are often referred to as 
Delaunay canonical variables L , G, H and 1 , g , h. The Keplerian elements are 
related to the Delaunay variables in the following way 


Pi 

= L 

(M3) 1 2 , 


= 1 

- M, 



p 2 

G 

= L (l ~ e 2 ) 1 2 , 

^2 

r g 

— CO , 

> l 

(30) 

P3 

= H 

= G cos i , 

^3 

= h 

= ft 

J 




Brouwer solved Hamilton’s equations using the disturbing function given in 
equation (26) by the method of canonical transformations. Among other things, 
he obtained the long-period and secular series expressions for the elements 1 , 
g and h as a function of the quantities J 2 , J 4 , a" , e" and i" etc. 2 The 
quantities a", e" and i" are the Brouwer mean elements. The significant 
terms of the secular disturbing function due to the earth, as derived by Brouw^er, 
can be found in Felsentreger's paper. He calls this function F s .3 This paper 
will call the secular part of the disturbing function due to the earth’s nonspherical 
gravitational field R s (E) . The quantity 

- J 2 , , 

R S(E) “ 4L 6 (l-e 2 ) 3 2 l 1_3cos2i ) 


+ 



_3 1 

128 (l -e 2 )5 2 


(5 “ 18 cos 2 i + 5 cos 4 i) 


X W. M. Kaula, ’’Celestial Geodesy,” NASA TN-D-1155, (March 1962), 8. 
^Brouwer, ’’Artificial Satellite Theory,” pp. 393-395. 

3 

Felsentreger, ’’Motion of Relay II,” p. 5. 
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_3_ 1 

32 (1 - e 2)« 2 


(l - 6 cos 2 i + 9 cos 4 i) 


15 1 

128 (i - e 2 ) 7 2 


(l - 2 cos 2 i - 7 cos 4 i) 


3 

128 



_5 3_ 

e 2 ) 7 2 " (1 -e 2 ) 



3 ~ 30 cos 2 


+ 35 cos 4 i) . 


(31) 


Another set of formulas which are very useful to describe the motion of a 
satellite is Lagrange's planetary equations. These equations give the time 
derivative of an orbital element as a function of the values of the six orbital 
elements. For the angular variables g and h, these equations are 


dg _ 

at " 


- cos i 

na 2 /l - e 2 sin i 


dR fl ~ e 2 dR 
di + na 2 e * e 


(32) 


and 


dh 

dt 


na 


2 VT 


TTT 


dJR 

di 


sin i 


(33) 


These equations can be simplified by defining a convenient set of canonical 
units of time and distance. Such a consistent set of units also facilitates all 
computations. By Kepler's third law, 

/jl - n 2 a 3 


(34) 
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where 


“ ^ (^e arth ^satellite) 

The quantity n is termed the mean angular motion of the satellite, where 

n z 2r Period . 

If one uses the following canonical units for time and distance, 

canonical time unit = 806.81242 seconds 


and 


then 


canonical length unit = 6378.166 kilometers, 


u n 2 a 3 1 . 

These units are currently used at Goddard Space Flight Center. With these units, 
the terms n 2 a 4 and na 2 in Lagrange’s planetary equations reduce to a and 


/a respectively. 


PART II 

DESCRIPTION OF THE MOTION IN THE ARGUMENT 
OF PERIGEE AND THE LONGITUDE OF THE ASCENDING NODE 
FOR THE EARTH SATELLITE RELAY II 
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CHAPTER I 

NATURE OF THE PROBLEM OF RESONANCE IN THE LUNI-SOLAR 
DISTURBING FUNCTION FOR RELAY II 

There is a computer program used at Goddard Space Flight Center which 
computes solar and lunar perturbations on a close earth satellite. The com- 
puter program is named SLOPE. It was formulated in a paper by Murphy and 
Felsentreger.l In general terms, the program employs the Lagrange planetary 
equations to form the derivatives of the orbital elements. The integration is 
then performed analytically to obtain the perturbation. For Relay II, the program 
predicts perturbations that are excessively large. Two terms in the luni- solar 
disturbing function can easily be identified as causing the problem. These 
terms are referred to as resonant. The resonant terms contain factors like 
cos 2 (g + h ). In the process of obtaining the argument of perigee, the factor 
is integrated and it yields the expression 

sin 2 (g +h~) 

2 (g +fi) 

^Murphy and Felsentreger, ’’Lunar and Solar Effects,” pp. 1-46. 
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The problem with this expression is that the time derivative of g and h 
have about the same magnitude but opposite sign. For Relay II, 

li = -1.1046373 degrees/day 
g = +1.1063814 degrees/day.* 

The term ( g + h ) is equivalent to In each analytical term where the time 
derivative of the longitude of perigee is the only term in the denominator, it 
produces a near singularity. 

The period for the longitude of perigee is excessively long, i.e., the 

period of 25 = 360/£ r 360/(0 +«) = 570 years . 

This means that the phenomenon of resonance occurs in the orbit of Relay II 
because of the derivative of the mean motion of an argument is equal to zero. When 
the derivative of the longitude of perigee 


''O 

Cjl) 


n + 


a: 


(35) 


then 


co - constant . 


This condition can be visualized by looking at figure 7. Physically, resonance 
occurs for Relay n (a close earth satellite, a = 1.7449 and a satellite of large 
eccentricity, e = .239) because it happens to be near a critical inclination. 

When the sun, moon and oblateness act simultaneously on such a satellite near a 


1 


Felsentreger, ’’Motion of Relay II,” 


pp. 26-23. 
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critical inclination, then resonance occurs. An appropriate value of the 
resonant inclination can be found by assuming that the main effects on g and 
h for Relay II are due to the oblateness effect, the J 2 or A 2 term. From 
Kozai's paper, 


and 



n cos 


i 



(36) 


(37) 


then 


• • 


n + 


A 2 / 5 . a . .) 

— 2 n ^2 ~~2 sin 1 2 i _ cos i J 


0 . 


Here p equals [a (l - e 2 )J . This reduces to the quadratic equation 


5 cos 2 i - 2 cos i ~ 1 0 


1 Y. Kozai, ’’The Earth's Gravitational Potential Derived from th^ Motion of 

Satellite 1958 Beta Two,” SAO Special Report, Number 22 (March, 1559). 2. 
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which has the solution 


cos i = - .2896, .6896 
or 

i = 106?833, 46?400 } 

The inclination of Relay II averages 46.32 degrees which is close to the latter 
value. 

m 

Since the orbit of Relay II is near this critical inclination, problems result 
when one analyzes the motion of the Keplerian elements for Relay II. Felsen- 
treger has investigated the elements e and i. This work will investigate the 
angular variables g and h. 


1 Peter Musen, private communication, 1968. 


CHAPTER H 

INPUTS FOR THE SOLUTION OF THE PROBLEM 

In order to obtain a solution describing the motion of g and h for the satel- 
lite Relay II, two inputs were required. One input was an appropriate analytical 
disturbing function R. The second input was a good stretch of observation data, 
which consisted of observed orbital elements of the satellite Relay II for a few 
hundred days. The disturbing function R is composed of several distinct com- 
ponents. Two components, R R(SM) and R g(E) , have already been mentioned. 

For R s , only the terms with coefficients J 2 , J 2 and J 4 will be considered. 

In addition to these components, there is a component R s (SM) which can be 
viewed as a secular term produced by the sun and moon. 

The expression for R s (SM) in Felsentreger’s paper is appropriate but a 
coordinate transform must be performed on the quantity i c .* The quantity i c 
cannot be considered constant over a long data arc. The quantities i" and i v 
are nearly constant and w." is almost linear with time. These three variables 
are more appropriate for problems which have long data arcs. From figure 8, 
the relationship between these spherical angles is 

^Felsentreger, '’Motion of Relay II,” pp. 2-5. 
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cos (180 ~i c ) _ cos i cos i * sin i sin i cos C 

Then 

sin 2 i c ^1 - cos 2 i cos 2 i 

• 2 cos i cos i sin i cos i. sin i 


(38) 


- sin 2 i sin 2 i . cos 2 0 I . 

The final expression for R S(SM) is 

R s f sm ) “ 32 [ n 2 m . “ -3 sin 2 i . ) 

n c 2 m c (-1^3 cos 2 i cos 2 i - 6 cos i cos i sin i " sin i _ cos 

* 3 sin 2 i sin 2 i cos 2 . jj ^L 4 (2 4 3e 2 ) (l - 3 cos 2 i)j . (39) 

The final form of the total disturbing function R for the system composed of 
Relay II, the earth, sun and the moon is 

^ ^R(SM) ^S(E) ^S(SM) (40) 

The disturbing function is composed of a series of terms containing the 
variables and parameters g,h, a, e, i, 0 , i?", i", and i . To obtain the 
necessary observed input values for a, e , i , g and h, elements were taken 
from Felsentreger’s paper. 1 He obtained these by processing minitrack obser- 
vation data with a Brouwer theory program. Such a program assumed the obser- 
vation can be explained by a theory using the zonal harmonic terms J 2 , J 3 , J 4 , 


1 Ibid., pp. 18-29. 


.i 
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J 5 and J2 in the earth’s gravitational potential. By this method one obtains 
values of the orbital elements g and h which still have the effect of the sun 
and moon in them. Solar radiation and long-period luni-solar perturbations were 
then computed using the SLOPE program (not including the long- period resonant 
terms) and subtracted from the Brouwer elements. This yields the so-called 
corrected elements e , i , g and h . These corrected elements still contain 

c c c c 

the effects of secular and resonant terms due to the sun and the moon. The use 
of these corrected elements eliminates the effects of the long-period luni-solar 
terms in the elements. The long-period luni-solar terms in the disturbing func- 
tion can then be eliminated. This makes it easy to spot periodic variation in 
the elements due to the odd zonal harmonics of the earth and to determine a 
value for J 3 . The computation for J 3 will be made as a check to the work. 

The corrected elements can still be viewed as observed orbital element data. 
These elements are given at roughly weekly intervals after the epoch date, Jan- 
uary 21, 1964, 21 hours and 41 minutes U.T., for some 650 days. The parameters 
ft" , ft*, i" and i , can be obtained for the desired date by polynomials given in 
Murphy and Felsentreger’s paper. * With the observed orbital elements data and 
the disturbing function, R (a, e c , i c , g c , h c , ft" , Cl. , i", i @ ), the problem of 
computing the perturbations in g and in h could be undertaken. 


^Murphy and Felsentreger, ’’Lunar and Solar Effects,” pp. 46-47. 



CHAPTER HI 

NUMERICAL METHOD OF SOLUTION 


Since the analytical series method as specified by the SLOPE computer pro- 
gram yielded excessively large perturbations for g and h , another approach 
was attempted. A better approach to computing the perturbations would be to 
avoid forming terms which had the derivative of the longitude of perigee in the 
denominator. This second approach to the problem again used Lagrange’s plane- 
tary equations to obtain the expressions for dg/dt and dh/dt. Upon evaluation 
of these derivatives at each observation date, it was found that the scatter in the 
magnitude of dg /dt was roughly twice the scatter in the amplitude of dh /dt . 

For example, from 87 to 101 days past epoch, dg/dt changed by 1.74 x 10~ 4 
degrees/day while dh/dt changed by .42 x 10 -4 degrees/day. Refer to the data 
in tables 1 and 2. These points are plotted on graphs 1 and 3. 

The integration of the derivatives dg/dt and dh/dt was done numerically 
using the trapezoidal rule for 84 observation dates. The formula used to obtain 
the perturbations in the argument of perigee was 


n-1 


Sg c 


ompu ted 





g (tj tl > + g (»!> 


<‘i.i - i i> 


(41) 
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where 

n =1,2,. . . ,84 

t 0 = start time of integration was 1/24/64, 21 hours, 41 minutes U.T. , 

and 

t i = i th observation date. 

Upon integrating dg/dt and dh /dt, the perturbations kg (t - t 0 )and bh(t 
were obtained at roughly weekly intervals for 655 days. 

The perturbations 8 g(t - t 0 ) increased almost linearly from 7.7 degrees 
at ten days past epoch to 720 degrees on day 654. The perturbation 8 h (t - t 0 ) 
decreased almost linearly by almost the same magnitude as S g. Refer to the 
data in tables 1 and 2. The computed values for the elements g and h were ob- 
tained by adding the computed perturbations Sg and Sh for any date past epoch 
to the value g c and h c on the date (t 0 ) at which the integration was begun. Then 

8ca«p„.ed <‘n) = 8 (»0> + <*„> 

and 

h . (t ) = h (t n ) + 8h (t ) 

computed v n / v 0 7 v n / 

where 

n = 1, 2, . . . , 84. 

To investigate the accuracy of the computed g and h, their values were 
compared to the input g c and h c values. The latter are essentially observed 
elements. Residuals were formed. Let 

^res ^ ^ n ^ — ^computed ” ^corrected 

and 


^*res ^ ^ n ^ ^computed ^ ^ n ^ ^corrected ^ n ^ 


(43) 
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The quantity g decreases almost linearly from zero to .34 degrees after 654 
days. During the same period h res increases almost linearly to .17 degrees. 

See tables 3 and 4 for the values of the compute, corrected and rt idual values 
of the elements. Plotting the values of the residuals against time for 654 days 
gives a nearly straight line. For g res , the line has noticeable arcs superimposed 
on it. For h res , the line is very straight and has only slight arcs on it. The 
quantities g rps and h ^ s are plotted on graphs 5 and 6. 

The main linear characteristic of the graphs for g, es and h res can be 
accounted for. Such secular deviation can arise because of errors in J 2 and 
J 4 . It also can arise because of omission of higher order even zonal harmonics 
such as J 6 , J 8 , etc. To better understand how the linear deviation of the com- 
puted from the observed elements arose, the residuals g and h were 
fitted to a least squares fit of the form 

g (t ) ~ A n + A. (t ) 

b res v n / 0 l N n / 

and 

h (t ) £ B ft + B. (t ). 

The resulting slope for g res had almost twice the absolute magnitude of the 
slope for h . By adjusting the zonal harmonic J 4 from -2.123 x 10 -6 to 
-1.672 x 10' 6 the magnitude of the slope for g res became almost the same as 
the slope for h res . Three values of J 4 were tested: -1.672 x 10' 6 , -1.840 x 
10~ 6 and -2.123 x 10" 6 . The least squares slope for g res in degrees per day 
for these J 4 values was -2.748 x 10~ 4 , -3.644 x 10 " 4 and -5.163 x 10~ 4 . The 
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corresponding slope for h res in degrees per day was 2.919 x 10 -4 , 2.781 x 10' 4 
and 2.489 x 10 ~ 4 . Such errors in J 4 are plausible. Higher order even har- 
monic.* presumably contribute to a lesser degree to the secular deviation. 

To further investigate the residuals between the computed and observed ele- 
ments, the straight-line least squares solution was subtracted from the residuals. 
This essentially left only the periodic components of g res and h res . Let these 
periodic components be denoted g prcs and h pres . See the data in tables 5 and 7. 
The equations for g and h are 

gpre, - 8.0, <‘„> * [A 0 + A 1 (V 1 (44) 

and 

h (t ) = h (t ) - [B n + B. (t )] . (45) 

Looking at the graphs of g pres and h s , graphs 7 and 10, which consist of 
84 points scattered from 10 to 654 days past epoch, one notices immediately 
that the graph for g pres is definitely periodic. The amplitude is about 27 seconds 
of arc or .0075 degrees. The maxima occur at about 319 and 633 days past 
epoch and minima occur at about 158 and 458 days past epoch. This indicates the 
presence of a function composed predominantly of a term which has a frequency 
of g (or a period of 338 days). 

The graph for h res appears to be quite different from the graph for g 

P res p res 

A long-period luni-so^ar term still seems to be present in h p . The luni-solar 
term has the same behavior as in Felsentreger's paper. 1 The amplitude of the 

^Felsentreger, "Motion of Relay II," p. 37. 
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term is about 8 seconds of arc. The existence of the term implies that the ele- 
ments e c , i c , h c , and g c , which were input to the program, were not corrected 
for this term. Since the amplitude is so small, a similar luni-solar term would 
only slightly effect the graph for g. The graph for h pres has an amplitude of 
only a few seconds of arc. It has many more arcs superimposed on it than g pres . 
This means that terms with a frequency 2h have a larger amplitude than terms 
with a frequency h . Also a great deal of noise is evident in this small amplitude 
data. To investigate the components present in g pres and h pres will require 
different approaches due to the different characteristics of the graphs. 

It is evident from the graph for g pres , that there is a fair amount of noise 
present in the data. Visual inspection of the curve indicates that it must have an 
amplitude less than 43 seconds or arc or .012 degrees. Four of the 84 points 
for g pres lie outside this limit. A typical maxima occurred around days 333 
and 347. The amplitude of g is .0079 and .0074 degrees on those date3. 

These values give a good indication of the actual amplitude of the periodic terms, 
g . The amplitudes of the periodic components in g had to be accurately 
determined to obtain a meaningful value for J 3 . This implied that the effect of 
the noise in the data for g had to be lessened. Either of two methods would 

pres 

be adequate: points which were fairly detached from the curve could be discarded 
or each pair of successive points could be averaged. The effect of both methods 
would be to lessen the maximum amplitude of g b . Since there was a great 
deal of noise, the first method was chosen. Four points which had an amplitude 
that was greater then .012 degrees and which would heavily bias a least squares 
solution were discarded. The remaining 80 points were fitted to an equation of 


the form 
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g 


p r <* s 


(t n ) - A 2 sin g (t n ) 


t A 3 cos g (t n ) 


* A 4 sin 2g (t n ) + A- cos 2g (t n ). (46) 


Then 


A 0 sin g 




COS £ 


(t n ) 


A^ i A? "sin 




Aj T A 3 cos 



+ 



The size of i A* + A 3 indicates the magnitude of the cos gt term which is pro- 
duced by the earth's J 3 zonal harmonic. This term was deliberately not in- 

Y~ 

eluded in the constructed input disturbing function. The magnitude > A 4 + A s 
indicates the size of the (J 4 / J 2 ) sin 2 g t term. Performing the least squares 
fit yields the following values for the least squares coefficients: 


A 2 - - .002649 degrees , 
A 3 - + .005335 degrees, 

A 4 - - .000717 degrees, 
A 5 = - .000532 de grees , 


i AJ + A 3 = .005956 degrees, 


v A 4 + Ag = .000893 degrees . 


Refer to graph 8 and table 6 which were constructed using these coefficients. 
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From the least squares value for A 1 2 +• A 2 and a formula used in Kozai's 
paper, a value for J 3 can be determined. 1 The perturbation in g due to the J 3 
harmonic by observation is 

' g, ( A; t A* cos g . (47) 

The perturbation according to Kozai is 

( - J 3 J 2 ) (1 - cos 2 i *■ o 2 cos 2 i) 

■ g | - cos g. (48) 

2 ca (1 - o 2 ) sin i 

The right-hand portions of each of the above equations are equated and solved for 
J 3 on each of the 80 observation dates. The average J 3 determined by using 
the angular variable g is -2.589 x 10 ' 6 . To get an idea as tc the error in this 
measurement, one can solve the case where no points were discarded, all 84 
points were used. Then 


A 2 r A 2 - .006508 degrees 

and 

• A 2 + A 2 - .001085 degrees . 

The quantity J 3 would then equal -2.826 x 10~ 6 . If the four g pres points are 
discarded as stray maxima points (noise), then the value of J 3 determined from 

1 Y. Kozai, ’’Second Order Solution of Artificial Satellite Theory without Air 

Drag,” Astronomical Journal, LXVII (September, 1962), 457. 


the angular variable agrees well with the J 3 determined from a study of the 
eccentricity by Felsentrcger. * His value for J was -2.503 x 10' 6 . 

Due to the presence of the conspicuous luni-solar term and of a great deal 
of noise, a determination of J 3 from b s posed difficulties. Also if one com- 
puted the predicted amplitude of the perturbation in h due to the J 3 harmonic 
term from the equation in Kozai’s paper, 


-b 


■>.i 


f “ J 3 J 2 ) «* ‘'os i 

2 a (1 - c 2 ) sin i 


(49) 


then the amplitude of the perturbation was on the order of a half second of arc. 2 
The magnitude of the scatter in the data for h s had about the same size. Thus 
a determination of J 3 from h^ res would have been meaningless. 


^Felsentreger, "Motion of Relay II " p. 9. 

“Kozai, "Second Order Satellite Theory," p. 457. 


CHAPTER IV 
CONCLUDING REMARKS 

The problem o* resonance complicating the computation of the angular or- 
bital elements for Relay II has been avoided. The problem had first been noticed 
when the computed values of the luni- solar perturbations for g and h using C e 
SLOPE method, as given by Murphy and Felsentreger, wore excessively large. 

In tV SLOPE method, the derivative of the angular variable was obtained by 
using the Lagrange p’anetary equations. These derivatives were then integrated 

analytically. The perturbation in g or h is then in the form of a series of 

, • 

trigonometric terms. Several of these terms have the factor ( g + h ) in the 
denominator. Since g and h have about the same magnitude but opposite sign, 
terms which contain their sum in the denominator are nearly singular. 

To avoid the problem encountered by the SLOPE method and to show the 
feasibility of computing the perturbations in the angular orbital elements 
accurately, another disturbing function of only a limited number of terms was 
formed. A coordinate transformation was done on i c to eliminate it from the 
disturbing function. Lagrange's equations were again formed to yield the de- 
rivatives of the angular elements. No term in the analytical expression for g or 
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h contained the factor (g + h ) in the denominator. The derivatives were evalu- 
ated and integrated numerically. The elements computed for g and h differed 
by approximately .3 degrees and .2 degrees from those observed at the end of a 
two-year period. The slope of the linear deviation in the computed and observed 
values for g and b can be added as a correction to any model attempting to ex- 
plain the motion of g and h for Relay II. After accounting for the secular portion 
of the difference between the observed and computed differences in the angular 
variables, the remaining periodic difference for g was shown to be mainly due 
to the J, zonal harmonic of the earth. The remaining periodic difference for h 
was mostly noise superimposed upon a small amplitude long period luni-solar 
term. The accuracy of this numerical approach to calculate the perturbations in 
g and h was tested by making a computation for J 3 using g. The value of J 3 
that was determined was consistent with the value obtained by Felsentreger, who 
studied the eccentricity of Relay II. 

The primary motive of this work was to show that the perturbations in the 
argument of perigee and the longitude of the ascending node for Relay II could be 
obtained by numerical integration of Lagrange's planetary equations. Since this 
method yielded an acceptable value for J 3 , the method was shown to be workable. 
The problem can now be reversed by solving only for the perturbations in g and 
h , rather than for the coefficient J 3 . To evaluate the perturbations, the analytical 
expression for the J 3 component of the disturbing function should be added to the 
constructed disturbing function R. Lagrange's planetary equations can then bo 
formed for the derivatives of the angular orbit elements. These derivatives 
should then be evaluated and integrated numerically. The resulting perturbations 
can then be used to generate an ephemeris for g and h. 
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a = semimajor axis of the satellite orbit 
e = eccentricity of the satellite orbit 
i = inclination of the satellite's orbital plane to the earth's 
equatorial plane 

M = mean anomaly of the satellite 

oj = argument of perigee of the satellite’s orbit 

ft = longitude of ascending node of the satellite's orbit 

f = true anomaly of the satellite 

n = mean angular motion of the satellite 

T = time of perigee passage 

£ = longitude of perigee 

t * • • 

a, e, i , co, ft, a>, l = time derivatives of the orbital elements 

8 (to) • h (t 0 ) = value of g ( or h) on January 24, 1964, 21 hours, 41 
minutes U.T. 

pj , P 2 . P 3 , , q 2 , q 3 = are the canonical variables. The are the canonical 

conjugates of the q i . 

P = orbital period 

•£. g, h = Delaunay canonical variables which are equivalent to 
the Keplerian elements M , o> and ft 
L = /a 
G = L /l - e 2 



H = G cos i 
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t n " t 0 = days since January 24, 1964, 21 hours, 41 minutes U.T. 
i h c , g c = e ^ emen ^ corrected for solar radiation pressure and longf- 
period luni-solar effects, excluding near-resonant terms 
y = direction of the vernal equinox 
m c = mass of sun (moon)/mass of sun (moon) + mass of earth 

n^, n c = mean motion of the sun (moon) relative to the earth 
i 0 , i c = inclination of the sun’s (moon's) orbital plane to the earth’s 
equatorial plane 

0 C = longitude of the mean ascending node of the sun’s (moon’s) orbit 
on the earth’s equator , measured from the mean equinox of date 
i " = inclination of the moon's orbital plane to the ecliptic = 

5.1453967 degrees 

0" = longitude of the mean ascending node of the lunar orbit on the 
ecliptic, measured from the mean equinox of date. 

J 2 . J 3 , J 4 = zonal harmonic coefficients in the earth’s gravitational potential 
A 2 = Kozai's equivalent to (3/2) J 2 
T. J. k = Cartesian unit vectors in the x , y , z directions respectively 
V = gradient vector operator. In Cartesian coordinates it is 

J _ J (9 — 

<9x 1 + <9y J + <9z k 


V 2 = Laplacian scalar differential operator V • V. in Cartesian 

coordinates it is 

<9 2 a 2 ( j 2 

<9x 2 Jy 2 dz 2 
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<) R (to , <1 R r? i 


partial derivative of the disturbing function R with respect to 
<■' (or i ), holding all the other elements constant 


>g. b = size of me perturbation in g (orb) over the time interval t - t 0 
V = potential as used in physics (r = ■ , V = 0), (r - 0, V = - » ) 

U = potential as used in astronomy U ~ * V . U is also called the 
forcing function, (r ~ * . U ’ 0). (r 0, U " ~ x ) 
k 2 = universal gravitational constant 
F = Hamiltonian oi a dynamical system 


T 


KE 


F = 
R = 



p _ 

‘'rc sm ) 
^S(E) = 

^S(SM) 


k 2 (nij + m 2 ). In canonical units ^ = 1 
kinetic energy of satellite 

force exerted on a body as given by Newton’s laws 
disturbing function. The contribution to the forcing function of a 
body which makes the body deviate from simple two body motion 
sum of the two near- resonant terms in the luni- solar disturbing 
function 

secular disturbing function 

the two resonant terms in the luni -solar disturbing function 
secular part of the disturbing function due to the earth’s spherical 
harmonics 

secular part of the disturbing function due to the sun and moon 
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m j , m 2 

r 



r ~ 


r = 


r = 


" l 


« • 



X. y. z = 


12 


r = 



-i - 


\ — 


mass of body one (two) as implied in Newton's second law 
radius of the earth 
initial vector position 
position vector 

first time derivative of the position vector 

second time derivative of the position vector 

position vector of body one with respect to the center of mass of 

the system 

second derivative of the vector 7 X 

Cartesian coordinate components of r 

first time derivatives of the Cartesian coordinates 

position vector of body two with respect to body one 

magnitude of the vector r 

magnitude of the vector A 

position vector of the perturbing body, the sun or moon, with respect 
to the earth 

geocentric colatitude of the satellite 
geocentric longitude of the satellite 

declination or the angular distance above or below the celestial 
equator at which the satellite appears in the sky 


R( r ) = function dependent on the coordinate r 
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TABLE 1 

COMPUTED DERIVATIVES AND PERTURB A HONS FOR THE ARGUMENT OF 

PERIGEE 


Obs. 

Number of 

Derivative 

Computed 

No. 

Days Since 

of the Arg. 

Pert, in the 


Epoch Date 

of Perigee 
(Degrees/Day) 

Arg. of Perigee 
(Degrees) 

2 

10*0000000 

1* 1056376 

7*7393016 

3 

17*0000000 

1*1059399 

19*4767064 

4 

24*0000000 

1*1096347 

23*2160706 

5 

31*0000000 

1*1039241 

30*9974200 

A 

30*0000000 

1*1036793 

30*6960536 

7 

49*0000000 

1*1063374 

46*4363096 

0 

92*0000000 

1*1056471 

54*1 796697 

9 

99*0000000 

1*1030746 

61*9191611 

10 

66*0000000 

1*1065963 

69*6544600 

11 

73*0600000 

1*1059902 

77*3930446 

12 

60*0(100000 

1*1067444 

69*1337123 

13 

67*0000000 

1*1067603 

92*8741606 

1A 

101*0000000 

1*1039700 

106*3539429 

13 

100*0000000 

1*1066337 

116*0933669 

* A 

119*0000000 

1*1037196 

123*6332026 

IT 

129*0000000 

1*1036936 

139*3133240 

10 

131*6000000 

1*1057796 

163*6397656 


1S6«6OOOOC0 

1*1057662 

171*3604160 

20 

165*0000000 

1*1030513 

173* 1 21 5220 

21 

172*0000000 

1*1096790 

166*0631 697 

22 

186*0000000 

1*1036399 

202*3455963 

23 

193*0000000 

1*1096092 

210*0065479 

24 

200*0000000 

l*109e912 

217*0277130 

23 

>07*0000060 

1*1037796 

225*5607256 

2 A 

214*0000000 

1*1066426 

233*3095703 

2T 

221 * 0000000 

1*1057756 

241*0504150 

2ft 

220*0000000 

1*1060369 

246*7912990 

29 

233*0000000 

1*1057662 

256*5319624 

30 

256*0000000 

l*106e044 

279*7543945 

31 

263*0000000 

1*1096674 

267*4953613 

32 

270*0000000 

1*1036111 

295*2365723 

33 

277*0000000 

1*1059160 

302*9777632 

34 

204*0000000 

1*1036273 

310*7192303 

3» 

291*0000000 

1*1036645 

310*4602051 

30 

296*0000000 

1*1037329 

326*2011719 

37 

305*0000000 

1*1097596 

333*941406? 

30 

312*0000000 

1*1062540 

341*6023730 

39 

319*0000000 

1*1097377 

349*4230997 

40 

326*0000000 

1*1092190 

357*1630164 

41 

333*0000000 

1*1056536 

364*9040527 

42 

340*0000000 

1*1036273 

372*6445312 


Epoch Date * 1/21/64, 21 hours, 41 minutes, U.T. 
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TABLE 1 — Continued 


> 


Obs. 

Number of 

Derivative 

Computed 

No. 

Days Since 

of the Arg. 

Pert, in the 


Epoch Date 

of Perigee 

Arg. of Perigee 



(Degrees/Day) 

(Degrees) 

43 

347*0000000 

1*1057196 

360*3650096 

44 

354*0000000 

1*1057577 

388* 1 254663 

45 

361*0000000 

1*1056661 

395*6657227 

46 

366*0000000 

1*1057567 

403*6059570 

47 

375.0000000 

1*1057166 

411*3461914 

4 « 

362 * noooooo 

1*1057720 

419*0866699 

44 

389*0000000 

1*1057663 

426*6271464 

SO 

395*0964365 

1*1057425 

433*5663594 

51 

402*0964355 

1*1056552 

441*30859.38 

52 

409*0964355 

1*1057606 

449*0466261 

53 

423*0964355 

1*1057976 

464*5302734 

54 

430*0964355 

1*1059437 

472*2717265 

55 

437*0964355 

1*1059017 

460*0134277 

56 

444*0964355 

1*1060905 

407*7553711 

57 

451*0964355 

1*1061144 

495*4982910 

56 

456*0964355 

1*1060953 

503*2412109 

50 

465*0964355 

1*1060162 

510*9636667 

60 

472*0964355 

1*1059046 

516*725e301 

61 

479*0964355 

1*1059006 

526*4672652 

62 

466*0964355 

1 « 10595 99 

534*2069644 

63 

493*0964355 

1*1060343 

541*9511719 

64 

500*0964355 

1*1059532 

549*6933594 

65 

507*0964355 

i <10600 es 

557*4353027 

66 

51 4.0964355 

1*1059675 

565*1774902 

67 

521*0964355 

2*1060104 

572*9196777 

66 

526* 0964355 

1*1060324 

580*6621094 

60 

535*0964365 

1*1060295 

586*4045410 

70 

549*0964355 

l*106OSei 

603*8694043 

71 

556*0964355 

1*1 061239 

611*6323242 

72 

563*0964355 

1 *1060065 

619*3750000 

73 

570*0964355 

1*1060467 

627*1 174316 

74 

577.0964355 

1*1060257 

634*6596633 

75 

504*0964355 

1*1060161 

642*6022949 

76 

591 .0964355 

1*1061344 

650*3449707 

77 

596*0964355 

1*1060696 

658*0876906 

74 

605*0964355 

1*1099626 

665*8300781 

70 

61 2*0964355 

1*1060429 

673*5720215 

60 

619*0964355 

1*1060646 

681*3146973 

61 

626*0964355 

1*1060591 

669*0576172 

62 

633*0964355 

1*1061192 

696*8005371 

63 

640*0964355 

I*106.07e2 

704*5434570 

84 

647*0964355 

1*1059570 

712*2851563 

65 

654*0964355 

1*1059665 

720*0273436 
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TABLE 





COMPETED DERIVATIVES 
THE LONGITUDE OF 

AND 

THE 

PERTURBATIONS 
ASCENDING NODE 

FOR 

Obs. 

Number of 


Derivative 

Computed 

No. 

Day: Since 
Epoch Date 


of the Long, 
of the Asc. Node 
(Degrees/Day) 

Pert, in the 
Long, of Asc. Node 
(Degrees) 

2 

10*0000000 


-1* 1040640 

-7*7285700 

3 

17*0000000 


— 1 • 1040298 

-19*4670990 

4 

24*0000000 


-1*1040773 

-23*1896669 

9 

31.0000000 


-1*1040268 

-30*9142609 

9 

30*0000000 


-1*1040850 

-38* 6426833 

7 

45*0000000 


-1*1040449 

-46*3715210 

A 

92*0000000 


-1*1040735 

-54*1001740 

9 

59*0000000 


-1*1040444 

-61*8287611 

10 

66*0000000 


-1*1040697 

-69*9573730 

11 

73.0000000 


-1*1040649 

-77* 266041 3 

12 

80*0000000 


-1*1041193 

-85*0148621 

13 

87*0000000 


-1*1041250 

-92*7438966 

14 

101.0000000 


-1*1040831 

-108* 2016280 

15 

108*6000000 


-1*1041126 

-115*9306940 

16 

1 15.0000000 


-1*1041431 

-123*6597900 

17 

129.0000000 


-1*1041489 

-139*1182550 

ie 

151* 0000000 


-1*1042042 

-163*4108429 

1* 

158*0000000 


-1*1042194 

-171*1405334 

2 0 

165*0000000 


-1*1042585 

-178*8703461 

21 

172*0000000 


-1*1042728 

-186*6004181 

22 

186*0000000 


-1*1042757 

-202*0606842 

23 

193*0000000 


-1*1042643 

-209*7907715 

2.4 

200*0000000 


-1*1043177 

-217*5210114 

29 

207.0000000 


-1*1042614 

-225*2512054 

26 

214*0000000 


-1*1043062 

-232*9813995 

27 

221.0000000 


-1*1042652 

-240*7115784 

29 

228*0000000 


-1*1043119 

-246*4418335 

29 

235.0000000 


-1*1042957 

-256*1721191 

30 

256*0000000 


-1*1043156 

-279*3630371 

31 

263*0000000 


-1*1043472 

-287*0937500 

32 

270*0000000 


-1* 1043434 

-294*6242166 

33 

277.0000000 


-1*1043758 

— *0S! a {*549316 

34 

284*0000000 


-1*1043539 

-310*2856445 

39 

291*0000000 


-1*1043625 

-318*0163574 

36 

298*0000000 


-1*1043348 

-325*7470703 

37 

305*0009000 


-1*1043482 

-333*4777632 

34 

312*0000000 


—1*1043844 

-341*2064961 

39 

319*0000000 


-1*1043501 

-346*9392090 

40 

326*0000000 


-1*1043596 

-356*6699219 

41 

333*0000000 


-1*1043177 

-364*4003906 

42 

340*0000000 


-1*1043653 

-372*1311035 
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TA BLK 2 — Continual 


Obs. 

Number of 

Derivative 

C -mputed 

No. 

Days Since 

of the Long. 

Pert in the 


Epoch Date 

of the Asc. Node 
(Degrees/Day) 

Long, of Asc. Node 
(Degrees) 

4* 

347.0600000 

-1.1043348 

-379 

8618164 

44 

354.0003000 

-1.1043367 

-387 

5922852 

45 

361.0000000 

-1.1043186 

-395 

3227839 

4* 

368.0000000 

-1.1043444 

-403 

0532227 

47 

375.0000000 

-1.1043472 

-410 

7839355 

48 

382.0000000 

-1.1043429 . 

-418 

5144043 

49 

3 0000000 

-1.1043520 

-426 

2451172 

50 

3"? 5.0964355 

-1.1043730 

-432 

9/77032 

51 

402.0964355 

-1.1043310 

-440 

7084961 

52 

409.0964355 

-1.1043587 

-448 

4392090 

53 

423.0964355 

-1.1043854 

-463 

9006348 

54 

430.0964355 

-1.1044445 

-471 

6310359 

35 

437.0964355 

-1.1044273 

-479 

3630371 

56 

444.0964355 

-1.1045084 

-4 87 

0944824 

57 

451.0964355 

-1.1045132 

-494 

8261719 

58 

458.0964355 

-1.1045256 

-502 

S561CS5 

59 

465.0964355 

-1.1045103 

-510 

2897949 

60 

472.0964355 

-1.1044903 

-218 

0217285 

61 

479.0964355 

-1.1044865 

-•25 

7531730 

62 

486.0964355 

-1.1045332 

-533 

4651074 

63 

493.0964355 

-1.1045475 

-541 

2170410 

64 

500.0964355 

-1.1045599 

-54e 

9469746 

65 

507.0964355 

-1.1045761 

-556 

6811523 

66 

51 4.0964355 

-1. 1045790 

-564 

4135742 

67 

521.0964355 

-1.1045856 

-572 

1457520 

68 

528.0964355 

-1.1046238 

-579 

8761738 

69 

535.0964355 

-1.1046124 

-587 

6105957 

70 

549.0964355 

-1.1046324 

•603 

0759277 

^1 

556. 0964355 

-1.1046696 

-61 C 

6oe593e 

72 

563. 0964355 

-1.1046410 

-618 

5412598 

73 

570.0964355 

-1.1046667 

-626 

2741699 

74 

577.0964355 

-1.1046486 

-634 

0068359 

75 

584.0964355 

-1.1046677 

-641 

7397462 

76 

591.0964355 

-1.1046839 

-649 

4726563 

77 

598.0964355 

-1.1046743 

-657 

2055664 

78 

605.0964355 

-1.1046362 

-664 

9364766 

79 

612.0964355 

-1.1046705 

-672 

6708984 

80 

619.096«*355 

-1.1046ei0 

-680 

4030066 

81 

626.0964355 

-1.1046734 

-688 

1369629 

82 

633.0964355 

-1.1046934 

-695 

6701172 

83 

640*0964355 

—1.1046877 

-703 

6030273 

84 

647.0964355 

-1.1046572 

-711 

3354492 

85 

654.0964355 

-1.1046562 

-719 

0663594 
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TABLE 3 

COMPUTED, OBSERVED AND RESIDUAL VALUES FOR THE 

ARGUMENT OF PERIGEE 



Number of 

Computed 

Observed 

Residual 


Days Since 

Argument of 

Argument 

(Comp.-Obs.) 


Epoch Date 

Perigee 

of Perigee 

(Degrees) 



(Degrees) 

(Degrees) 


to 

ooooooo 

198*7867444 

195*7661274 

-0*0113831 

17 

0000000 

209*4960632 

263*8079803 

-0*0119171 

24 

ooooooo 

211 *2954431 

211*2866376 

-0*0211948 

31 

ooooooo 

219*9747929 

216*9953461 

-0*0205536 

9a 

ooooooo 

2 26*7142181 

226* 7436965 

-0*0296783 

45 

ooooooo 

294*4536749 

234*4820099 

-0*0283356 

52 

ooooooo 

242*1930542 

242*2241669 

-0*0311127 

59 

OOOOOOO 

249*9925256 

249*9664917 

-0*0339661 

68 

ooooooo 

257*6716309 

267*7160645 

—0*0444336 

79 

OOOOOOO 

265*4111928 

265*4548340 

-0*0437012 

80 

OOOOOOO 

273*1506789 

273*1969746 

-0*0480957 

87 

OOOOOOO 

280*8913574 

280*9433594 

-0*0520020 

101 

OOOOOOO 

296*3710937 

296*4265137 

-0*0554199 

109 

ooooooo 

304*1105957 

304*1713867 

-0*0607910 

115 

ooooooo 

3ii *esose59 

311*9147949 

-0*0642090 

129 

ooooooo 

327*3305664 

32 7*4052734 

-0*0747070 

151 

ooooooo 

351*6569824 

351*7507324 

-0*0937500 

158 

ooooooo 

359*3977051 

359*4963379 

-0*0986328 

165 

ooooooo 

7*1389160 

7*2409e6e 

— 0*10207^ o 

172 

ooooooo 

14*3803711 

14*9842663 

—0* 10391 i. * 

196 

ooooooo 

30*3627930 

30*4726651 

-0*1100922 

193 

0000006 

38*1037598 

38*2146436 

-0*1110840 

200 

ooooooo 

45*8449707 

45*9591675 

-0*1141968 

207 

ooooooo 

53*5859375 

53*69e2eao 

-0*1123505 

214 

ooooooo 

61*3269043 

61*4506209 

-0*1239166 

221 

ooooooo 

69*0676270 

69*1942291 

-0*1266022 

228 

ooooooo 

76*8065936 

76*9362163 

-0*1276245 

235 

ooooooo 

84*5493164 

84*6751556 

-0*1258392 

256 

ooooooo 

107*7717285 

107*9191132 

-0*1473846 

263 

0900000 

115*5126953 

115*6349835 

-0*1422882 

270 

ooooooo 

123*2539062 

123*3968600 

-0*1449738 

277 

ooooooo 

130*9951172 

131*1443176 

-0*1492004 

284 

ooooooo 

1 38*7365723 

13 6 *6918762 

-0*1553040 

291 

ooooooo 

146*4775391 

146* 6333771 

-0*1556380 

298 

ooooooo 

154*2185659 

154*3786011 

-0*1600952 

305 

ooooooo 

161*9587402 

162*1191101 

-0*1603699 

312 

ooooooo 

169*6997070 

169*6631744 

-0*1634674 

319 

ooooooo 

177*4404297 

177*6011353 

-0*1607056 

326 

ooooooo 

185*1811523 

168*3511047 

-0*1699524 

333 

ooooooo 

192*9213867 

193*0941315 

-0*1727448 

340 

ooooooo 

200*6616852 

200*8407593 

-0*1788940 
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TABLE 3 — Continued 


Number of 
Days Since 
Epoch Date 


000000 
j **. 0000000 

16 1. 0000000 

169.0000000 
379*0000000 

392.0000000 
399*0000000 
395*0964359 
403*0964355 
409*0464355 
423*0964355 
430*0464355 
437*0964355 
444*0464355 
491*0964355 
458*0964355 
465*0964355 

472.0964355 
479*0964355 
466*0964355 
493*0964355 
500*0064355 

507.0964355 
514*0964355 
521*0964355 
526*0064355 
535*0964355 
549*0964355 

556.0964355 
563*0964355 

570.0964355 

577.0964355 

584.0964355 
591*0964355 
596*0964355 
605*0964355 
612*0964355 
619*0964355 

626.0964355 
633*0964355 
640*0964355 

647.0964355 
654*0964355 


Computed 
Argument of 
Perigee 
(Degrees) 

206*402343? 
216*1426223 
223* 6R3C566 
231*6232910 
239*3 635254 
247*1040039 

254*8444524 
261*8656934 
269*325927? 
277*0661621 
292*5476074 
300*2890625 
306*0307617 
315*7727051 
323*5156250 
331 *2565449 
339*0012207 
346*7431641 
354*4846101 
2*2263184 
9*9685059 
17*7106934 
25*4526367 
33*1946242 
40*9370117 
48*6794434 
56*4218750 
71*9067383 
79*6496582 
07*3923340 
95*1347686 
102*8771973 
110*6196289 
118*3623047 
126*1052246 
133*8474321 
141*5693555 
$49* 332031 2 
157*0749512 
164*8176711 
172*5607910 
180*3924902 
186*0446777 


Observed 
Argument 
of Perigee 
(Degrees) 

206*5828094 
216*3298593 
224*0752106 
231*8216246 
239*8066609 
247*314986? 
285*0524717 
201*8024902 
269*5471191 
277*2915039 
292*7868652 
300*8310659 
306*2761230 
316*0175781 
323*7636719 
331 • 51 46484 
334*2566152 
347*0026358 
354*7470703 
2*4696586 
10*2461361 
17*9869221 
25*7221647 
33*4796889 
4^*2181702 
46*9666737 
56*7064198 
72*1979210 
79*9532013 
17*6879272 
95*4409027 
103*1825104 
110*9314423 
119*6692200 
126*4139453 
134*1591034 
141*9064454 
149*6846920 
157*3941040 
165*1446078 
171*6916931 
180*6322021 
168*3654623 


Residuel 

(Comp.-Obs.) 

(Degrees) 

-0*1604657 
-0*1667371 
-0*1921839 
-0*190333? 
-0*2031658 
-0*2109628 
—0*21 39b93 
-0*2167969 
-0*2211914 
-0*2253416 
-0*2392870 
-0*2419434 
-0*2483613 
-0*2440730 
-0*2460469 
-0*2861038 
-0*2543949 
-0*2895219 
-0*2624612 
-0*2635403 
-0*2716303 
-0*2782206 
-0*2755260 
-0*2650647 
-0*2611864 
-0*2892303 
—0. 2668446 
-0*2911635 
-0*3038431 
-0*2955933 
-0*3061371 
-0*3053131 
-3*31 16134 
-0*3069153 
-0*3081207 
-0*3116913 
-0*3170929 
-€•3226607 
-0*3191526 
-0*3267365 
-0*3309021 
-0*3297119 
-0*3407745 
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TABLE 4 

COMPUTED, OBSERVED AND RESIDUAL VALUES FOR THE 
LONGITUDE OF THE ASCENDING NODE 


Number of 

Computed 

Observed 

Residual 

Days Since 

Longitude 

Long, of 

(Comp.-Obs.) 

Epoch Date 

of Asc. Node 

Asc. Node 

(Degrees) 


(Degrees) 

(Degrees) 


10*0000000 

212*5613096 

212*5596907 

0*0014191 

L 7. 0000000 

204*6327769 

204*8291626 

0*0036163 

24*0000000 

197*1042175 

197*0984192 

0*0097963 

31*0000000 

169*3756256 

169*3675642 

0*0060414 

36*0000000 

161*6470032 

161*6370544 

0*0099467 

45*0000000 

173*9163655 

173*9066146 

0*0119509 

52.0000000 

166*1897125 

166*1775513 

0*0121613 

50*0000000 

158*4611053 

156*4484556 

0*0126499 

66*0000000 

150*7325134 

150*7146695 

0*0176239 

73*0000000 

143*0038452 

142*9832153 

0*0206299 

60*0000000 

135*2750244 

135*2530365 

0*0219679 

67*0000000 

127*5459900 

127*5232649 

i? *0227091 

101*0000000 

112*0860585 

112*0612640 

0 0267944 

106*0000000 

104*3591919 

104*3300323 

0 *0 291596 

115.0000000 

96*6300964 

96*6002608 

0 * 029 eiS 7 

124.0000000 

6 1*1716309 

81*1366272 

0*0350037 

151*0000000 

56*6790436 

56*8399200 

0*0391235 

156.0000000 

49*1493530 

49*1082153 

0*0411377 

165.0000000 

41*4195404 

41*3734569 

0*0460819 

172. 0000000 

33*6894664 

33*6444244 

0*0450439 

166.0000000 

16*2292023 

16*1796401 

0*0493622 

193.0000000 

10*4991179 

10*4461279 

0*0529900 

200 * 0000000 

2*7668760 

2*7165364 

0*0523396 

20 7.0000000 

35 S * 03 e 5742 

354*9651074 

0*0534668 

214.0000000 

347*3063496 

347*2529297 

0*0554199 

221.0000000 

339*5761250 

339*5207520 

0*0573730 

220*0000000 

331*6479004 

331*7676416 

0*0600566 

235*0000000 

324*1176758 

324*0571269 

0*0605469 

256*0000000 

300*9267576 

300*6618164 

0*0649414 

263*0000000 

293* 1960449 

293*1286621 

0*0673628 

270.0000000 

285*4655762 

265 • 3967285 

0*0666477 

277.0000000 

277*7340633 

277*6638164 

0*0710449 

264.0000000 

270*0041504 

269*9328613 

0*0712691 

291*0000000 

262*2734378 

262*1989746 

0*0744629 

296*0000000 

254*5426314 

254*4674225 

0*0704069 

305*0000000 

246*6121185 

246*7346039 

0*0775146 

312*0000000 

239*0814056 

229*0040894 

0*0773163 

319*0000000 

231*3506927 

231*2694702 

0*0812229 

326*0000000 

223*6199799 

223*5367126 

0*0832672 

333*0000000 

215*8895111 

215*8036346 

0*0656768 

340*0000000 

206*1567962 

208*0719604 

0*0666376 
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TABLE 4 — Continue' 


Number of 

Computed 

Observed 

Residual 

Days Since 

Longitude 

Long, of 

(Ccmp.-Obs.) 

Epoch Date 

of Asc. Node 

Asc. Node 

(Degrees) 


(Degrees) 

(Degrees) 

34 7*0000000 

200*4280853 

200*3385010 

0*0695644 

394*0000000 

192*6976166 

192*6063365 

0*0912781 

361*0000000 

184*9671478 

184*6746932 

0 *0922546 

368*0000000 

177*2366791 

177*1418610 

0*0946161 

373*0000000 

169*5059662 

169*4085063 

0*0974579 

362*0000000 

161*7754974 

161*6759644 

0*0995331 

369*0000000 

154*0447845 

153*9466553 

0*0961293 

395* 09643 53 

147*3121165 

147*2116663 

0*1002592 

402*09*4363 

t 39*5814056 

139*4793396 

0*1020660 

409*0964353 

131*8906927 

131*7459717 

0*1047211 

423*0964393 

114*3892670 

116*2817230 

0*1075439 

430*0964369 

108*65606^8 

106*5491*60 

0*1069476 

437*0964359 

100*9270020 

100*9151093 

0*1118927 

444*0964355 

93*1995566 

93*0623517 

0*1132050 

491*0964355 

85*4638672 

65*3509064 

0*1129606 

498*0964355 

77.7319336 

77*6152496 

0*1 166640 

465*0464355 

70*0002441 

69*8645673 

0*1156769 

472*0964355 

62*2663105 

62 *1 520644 

0*1162262 

479*0«64335 

54*5368652 

54*4187164 

0*1161486 

466*0964335 

46*8049316 

46*6655164 

0*1194153 

493*0964355 

39*0729980 

36*9510603 

0*1219177 

500*0964355 

31*3410645 

31*2163066 

0*1247559 

507*0564355 

23*e0eee67 

23*4606519 

0*1262349 

514*0464355 

15*8764648 

15*7472266 

0*1292362 

521*0964355 

8*1442871 

8*0136070 

0*1306801 

528*0464365 

0*4118652 

0*2791299 

0*1327353 

535*0964355 

352*6794434 

352*5441695 

0*1352539 

545*0564336 

337*2141113 

337*0761719 

0*1379395 

556*0564355 

329*4814453 

329*3403329 

0*1411133 

563*0564355 

321 *7467793 

321*6061543 

0*1406250 

570*0964355 

314*0156691 

313*8742676 

0*1416016 

S77.0464'55 

306*2832031 

306*1394043 

0*1437986 

584*0964355 

298*5502930 

298*4033203 

0*1469727 

551 *0964365 

290*8173828 

290*6701660 

0*1472166 

398*0964355 

283*0844727 

282*9362793 

0*1461934 

605*0964365 

275*3515625 

275*2001953 

0*1513672 

612*0944355 

267*6191406 

267*4675293 

0*1516113 

619*0464355 

259*8862305 

2! 9*7302246 

0*1560059 

626*0444365 

252*1530762 

251*9956512 

0*1574249 

633*0944355 

244*4199219 

244*2615356 

€•1563562 

640*0964355 

236*667011 7 

236*5265350 

0*1604767 

647*0964355 

228*9545898 

226*7689; 39 

0*1666169 

654*0964359 

221*2216797 

221*0537567 

0*1679230 
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TABLE 5 


RESIDUAL OF THE 

COMPUTED AND OBSERVED 

VALUE OF THE 

ARGUMENT OF 

PERIGEE AND ITS PERIODIC 

COMPONENT 

Number 

Residual of the 

Periodic 

of Days 

Comp, and Obs Value 

Component of 

Past Epoch 

of the Arg. of Per. 

Original Res. 


(Degrees) 

(Degree*) 

19.0000000 

-0*0113631 

0*0025341 

1 7.0000000 

-0*0119171 

0*0056146 

24.0000000 

-0*0211945 

-0*0000461 

31.0000000 

-0*0205536 

0*0042074 

38.0000000 

-0*0296763 

-0*0013026 

49.000^000 

-0*0283366 

0*0036546 

52.0000000 

-0*0311127 

0*0044921 

59.0000000 

-0*0339661 

0*0052533 

66.0000000 

-0*0444336 

-0*0015996 

73.9000000 

-0*0437012 

0*0027474 

eo.ooooooo 

-0*0460957 

0*0019675 

87.0900000 

-0*0320420 

0*0016759 

101.0000000 

-0*0334199 

. 0*0054671 

109.0000000 

-0*0607910 

0*0037306 

1 15.0000000 

-0*0642090 

0*0039272 

t 29.0000000 

-0.0747070 

0*0006584 

151.0000000 

-0*0937500 

-0*0070244 

X 58.0000000 

-0*0966328 

-0*0082926 

1 65. 0000000 

-0*1020708 

-0*0081160 

172.0000000 

-0*1039152 

-0*0063456 

1 96.0000000 

-0*1100922 

-0*0052935 

193.0000000 

-0*1110840 

-0*0026708 

2O0.6COOOO9 

-0*1141968 

-0*0021690 

207.0000000 

-0*1123505 

0*0032920 

214 .0000000 

-0*1239166 

-0*0046596 

221.0000000 

- 0*1266022 

-0*0037305 

228.0000000 

-0*1276245 

-0*0011383 

235.0900090 

-9*1250392 

0*0042616 

2S6.00OO900 

-0*1473846 

-0 *0064399 

263.0000900 

-0*1422862 

0*0022711 

270.0000000 

-0*1449736 

0*0032002 

277. 0 000000 

-0*1492004 

0*0025881 

204.0000009 

-0 • 1553340 

0*0000992 

291.0000000 

—0*1 358360 

0*0031797 

298.0000000 

-0*1600952 

0*0025371 

705.9000000 

-0*1603699 

0*0058771 

31 2.0000400 

-0*1614674 

0*0063941 

719.0000000 

-0*1607056 

0*0127706 

326.0000000 

-0*1699524 

0*0071384 

337.009009C 

-0*1727446 

0*0079606 

340*4000000 

-0*1766940 

0*0054259 


( 
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TABLE 5 — Continued 


Number 

Residual of the 

Periodic 

of Davs 

Comp, and Obs. Value 

Component of 

Past Epoch 

of the Arg. of Per. 

Original Res. 


(Degrees) 

(Degrees) 

3 % 7*0000000 

-€•1604657 

0*0074669 

354.0000000 

—0*1 667371 

0*0048121 

361.0000000 

-0*1 921539 

0*0030099 

36e.0000000 

-0*1963337 

0*0004447 

3 75.0000000 

-0*2031565 

-0*0007625 

362.0000000 

-0*2109526 

-0*0049452 

309.0000000 

-0*2139593 

-0*0043671 

395.0964355 

-0*2167969 

-0*0040267 

407.0964355 

-0*2211914 

-0*0046066 

409.09643*5 

-0*22 53 415 

-0*0053424 

423.0964355 

—0 *23 *12576 

-0*0120292 

430*0964355 

-0*2419434 

-0*0111001 

437.096^355 

-0*2453613 

-0*0109035 

444.0964^55 

-6*2446730 

-0*0068006 

451.0964355 

-0*2460469 

-0*0063596 

450.0964355 

-0*2561035 

-0*0106019 

465.0964355 

"0*2543945 

-0*0054763 

472.0964355 

-0*2595215 

-0*0069906 

479.0964355 

- 0*2624512 

-0*0063057 

466.0964355 

•>0.2635403 

-0*0037802 

493.0964355 

-0*2716303 

-0*0062556 

500.0964355 

-0c27522e6 

-0*0112395 

507.0964355 

-0*2755260 

-0*0049240 

514.0964355 

-0*2650647 

-0*0108461 

521 • 0964355 

-0* 2611554 

-0*0033253 

526.0964355 

-0*2692303 

-0*0077826 

535.0964355 

-0*2665445 

-0*0014624 

549. 0964355 

-0* 2911 835 

0*0011061 

556.0964355 

-0*3035431 

-0*0076369 

563.0964355 

-0*2955933 

0*0039275 

570.0964355 

-0*3061371 

-0*0030017 

577.0964355 

-0*3053131 

0*0014369 

564.0964355 

-0*3118134 

—0*0014466 

591. 0964355 

-0*30691 Sol 

0*0070639 

596.0964355 

-0*3061207 

0*0094731 

605.0964355 

-0*3116913 

0*0095171 

61 2.0964355 

-0*3170929 

0*0077301 

6: 9.0964355 

-0*3226607 

0*0055769 

626.0964355 

-0*3191528 

0*0128994 

633.0964355 

-0*3267365 

0.0C69304 

640.0964355 

— 0*3 30 9071 

0*0083793 

647.0964355 

-0*3297119 

0*0131841 

654.0964355 

-0*3407745 

0*0057361 
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TABLE 6 

PERIODIC RESIDUAL, LEAST SQUARES APPROXIMATION 
AND REMAINING RESIDUAL, FOR (g) 


Obs. 

Number of 

Periodic 

Least Squares 

Remaining 

No. 

Days Since 

Residual 

Value 

High Frequency 


Epoch Date 

(Degrees) 

(Degrees) 

Periodic Res. 
(Degrees) 

2 

19*0000000 

0*0025341 

0*0039983 

-0*0014642 

3 

17*0000000 

0*0056146 

0*0033963 

0*0022184 

4 

24*0000000 

*0*0000491 

0*0027798 

-0*0028279 

5 

31*0000009 

0*0092074 

0*0021614 

0*0020459 

6 

38*0000000 

*0*0013928 

0*0015504 

—0*0028532 

7 

45*0000000 

0*0036546 

0*0009532 

0*0027014 

8 

52*0000000 

0*0044921 

0*0003731 

0*0041190 

9 

59*0000000 

0*0052533 

-0*0001890 

0*0054423 

0» 

66*0000000 

-0*0015996 

-0*0007340 

-0*0008656 

a 

73*0000000 

0*0027474 

-0*0012643 

0*0040117 

12 

80*0900000 

0*0019675 

-0*0017825 

0*0037500 

13 

87*0000000 

0*0016759 

—0*0022909 

0*0039667 

14 

101*0000900 

0*0054e71 

-0*0032805 

0*0067676 

15 

108*0000000 

0*0037306 

-0*0037578 

0*0074684 

15 

115*0000000 

0*0039272 

-0*0042166 

0*0081438 

1 7 

129*0009000 

0*0006564 

-0*0050419 

0*0057003 

18 

1 51 *0900000 

*0*0070244 

—0* 0059780 

-0 *0011464 

19 

158*0009000 

*0*0082926 

-0*0059639 

-0*0023267 

29 

165*0000000 

-0*0081160 

—0*0059393 

-0*0021767 

21 

1 72*0000000 

-0*0063458 

-0*0057938 

-0*0005521 

22 

ie6* 0009000 

-0*0052935 

^0*0051154 

-0*0001761 

23 

193*0000000 

-C*0C2670e 

-0*004Sei7 

0*0019109 

24 

200*0000000 

-0*0021690 

-0*9039257 

0*0017568 

2 S 

297*0000000 

0*0032920 

-0*0031597 

0*0064516 

26 

214*0000000 

-0*0046596 

-0*0023002 

-0*0023593 

27 

221*0000000 

-0*0037305 

-0*0013685 

-0*0023620 

28 

228*0000000 

-0*00113e3 

-0*0003891 

-0*0007491 

29 

235*0000000 

0*0042616 

0*0006110 

0*0036506 

30 

256*0000000 

-0*9054399 

0*0034535 

-0*0098934 

31 

263*0000000 

0*0022711 

0*0042590 

-0*0019679 

32 

270*0000000 

0*0032002 

0*0049554 

-0*0017552 

33 

277*0000000 

0*0025861 

0*0055257 

-0*0029376 

34 

284*0000000 

0*0000992 

e*oo59seo 

—0 *0056588 

35 

291*0000000 

0*0031797 

0*0062457 

-0*0030659 

36 

298*0000000 

0*0025371 

0*0063875 

-0*0036503 

37 

345*0000000 

0*0C5e771 

0*0063874 

-0*0005103 

38 

312*0900000 

0*9063941 

0*0062541 

0*0001400 

39 

326*0000000 

0*0071384 

0*0056415 

0*0014969 

90 

333*0900000 

0*0079606 

0*0051955 

0*0027651 

41 

340*0000009 

0*0054259 

0*0046809 

0*0007451 


r 
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TABLE 6 — Continued 


Obs. 

Number of 

Periodic 

Least Squares 

Remaining 

No. 

Days Since 

Residual 

Value 

High Frequency 


Epoch Date 

(Degrees) 

(Degrees) 

Periodic Res. 
(Degrees) 

42 

347*0000000 

0*0074685 

0*0041161 

0*0033527 

43 

354*0000000 

0*0048121 

0*0035ie7 

0*0012934 

44 

361*0000000 

0*0030099 

0*0029041 

0*0001057 

45 

368* 0000000 

0*0004447 

0*0022853 

-0*0018406 

45 

375* 0009000 

-0*0007625 

0*0016722 

-0*0024347 

47 

382*0000000 

-0*0049452 

0*0010718 

-0*0060169 

45 

389*0000006 

-0*0043671 

0*00048ei 

-0*0048551 

49 

395*0964355 

-0*0040267 

-0*0000055 

-0*0040211 

50 

402*0964355 

-0*0046066 

-0*0005560 

—0 *0042506 

5! 

409*0964355 

-0*0053424 

-0*0010908 

-0*0042515 

52 

430*0964355 

-0*0111001 

-0*0026267 

-0*0064734 

53 

437* 0964355 

-0*0109039 

-0*0031201 

-0*0077834 

54 

444*0564355 

-0*006e006 

-0*0036021 

-0*0031965 

55 

451 «0964355 

—0*00 63556 

-0*0040678 

-0*0022920 

56 

458*0964355 

-o*cioeoi5 

-0*0045055 

-0*0062523 

57 

465*0964355 

-0*0054783 

-0*0049171 

-0 *0005612 

58 

472*0564355 

-0*0069906 

-0*0052780 

-0*0017126 

59 

479*0564355 

-0*0063057 

-0*0055781 

-0*0007275 

60 

486*056^355 

-0*0C37e02 

-0*0056026 

0*0020224 

61 

493*0564355 

-0*0082556 

-0*0059363 

-0*0023193 

62 

£00 •'09 643 5 5 

—0*0 1 1 2355 

-0*0059655 

-0*0052740 

63 

507* 0564355 

— 0* 1049240 

-0*00S87e2 

0*0009542 

64 

514*0564355 

-0*0106461 

-0*0056657 

-0*0051804 

65 

521*0504355 

-0*0033253 

-0*0053230 

0*0019977 

66 

528*0564355 

-0*0077826 

-0*0046497 

-0*0025326 

67 

535*0564355 

-0*0014824 

—0*0042504 

0*0027680 

66 

549*0564355 

0*0011061 

-0*0027169 

0*0038250 

69 

556*0964355 

-0*0076369 

-0*0018165 

-0*0058204 

70 

563*0564355 

0*0035275 

-0*0008564 

0*0047839 

71 

570*0964355 

-0*0030017 

0*0001374 

-0*0031391 

72 

577*0564355 

0*0014369 

0*0011371 

0*0002996 

73 

584*0964355 

-0*00144 6e 

0*0021141 

-0*0035629 

74 

591*0964355 

0*0070639 

0*9030406 

0*0040233 

75 

558*0564355 

0*0094731 

0*0038906 

0*0055625 

76 

605*0964355 

0*0095171 

0*0946410 

0*0048761 

77 

612*0964355 

0*0077301 

0*0052727 

0*0024574 

78 

619*0564355 

0*0055769 

0*0057714 

-0*0001945 

79 

633*0964355 

0*0069304 

0*0063385 

0*0025919 

80 

640*0964355 

0*0083793 

0*0064047 

0*0019747 

81 

654*0964355 

0*0057361 

0*0061344 

-0*0003983 
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TABLE 7 

RESIDUAL OF THE COMPUTED AND OBSERVED VALUE OF THE LONGITUDE 
OF THE ASCENDING NODE AND ITS PERIODIC COMPONENT 


Number 

Residual of the 

Periodic 

of Days 

Conip. and Obs. Value 

Component of 

Past Epoch 

of the Long, of Asc. N. 

Original Res. 


(Degrees) 

(Degrees) 

10*0000000 

0*0014191 

-0.0025982 

17.0000000 

0.0036163 

-0.0021439 

24.0000000 

C.0C57963 

—0.0017049 

31.0000*00 

C.0C6O414 

-0.0012048 

38.0000000 

C .0099467 

-0.0010404 

no ooooo 

0.0115509 

-0.0011811 

52.0000000 

0.0121613 

-0.0023137 

59. 0000000 

0.0126495 

-0.0035664 

66.0000000 

0.0176239 

-0.0003370 

72.0000000 

0.0206299 

0.0009260 

60.0000000 

0.0219879 

0.0005411 

67.0000000 

0.0227051 

-0.0004647 

101.0000000 

0.0267944 

0.0001188 

1 06.0000000 

0.0291595 

0.0C07409 

115.0000000 

€.0296157 

-0.0003459 

129.0000000 

0*0350037 

0.0013562 

151.0000000 

0.0391235 

-0.0000018 

i5e.ooooooo 

0.0411377 

0.0002695 

165.0000000 

0.0460eiS 

0.0034704 

1 72.0000000 

0 .04 50439 

0.0006898 

166.0000000 

0.0493622 

0.0015221 

193.0000000 

0.0529900 

0.0034070 

200.0000000 

0.0523396 

0.0010136 

207.0000000 

0.0534666 

0.0003979 

214.0000000 

0.C554199 

0.0006081 

221.0000000 

0.0573730 

0.0008183 

228.0000000 

0.0600566 

0.0017609 

235.0000000 

0.0 6054 69 

0.0005062 

256.0000000 

0.0649414 

-0.0003281 

263.0000000 

0.C673626 

0.0003704 

270.0000000 

©•06*8477 

0.0000923 

277.0000000 

0.0710449 

0.0005466 

2e4.0000000 

0.0712891 

-0.0009522 

291.0000000 

0.0 74 4629 

0.0004786 

298.0000000 

0.07540,69 

-0.0003183 

305.0000000 

0.0775146 

0.0000445 

312.0000000 

0.0773163 

-0.0018968 

319.0000000 

0.0612225 

0.0002665 

326.0000000 

0.0632672 

0.0005682 

333.0000000 

0.0656765 

0.0014345 

340.0000000 

0.0668376 

0.0006528 
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TABLE 7 — Continued 


Number 

Residual of the 

Periodic 

of Days 

Comp, and Obs. Value 

Component of 

Past Epoch 

of the Long, of Asc. N. 

Original Res. 


(Degrees) 

(Degrees) 

3*7* OQOCOOO 

0*0895844 

0*0016565 

354.ooooono 

0*0912781 

0*0016073 

361*0900000 

0*0922546 

0*0008409 

36l?.COP0O0<) 

0*0948181 

0*0016614 

375.0000PCP 

0*0974579 

0*0025582 

282*0000000 

0*09$5331 

0*0028905 

389*0000000 

0*0961293 

•0*0002563 

395*09643*5 

0*1002502 

0*0003467 

402*09643*5 

0*1020660 

0*0004196 

404.09643*5 

0*1047211 

0*0013316 

423.0964355 

0*1075439 

0*0006686 

430*09643*5 

0*1089478 

0*0003299 

437.0464355 

0*1118927 

0*0015316 

444* 0^64355 

0*1132050 

0*0011006 

431.0464355 

0*1129608 

-0*0008863 

455*0464355 

0*1166840 

0*0010939 

465. 0964355 

0*1156769 

-0*0016561 

4 7 2. 0 9643 £5 

0*1162262 

-0*0028498 

479. 0464355 

0*U6146e 

-0*0026701 

466*0964355 

0* 1194153 

-0*0031466 

403.0964355 

0*1219177 

-0*0023871 

SCO. 0964355 

0*1247559 

-0*0012919 

507.0964355 

0*1262349 

0*0004442 

514*09643*5 

0*1292362 

-0*0002974 

521.0964355 

041 30 6801 

-0*0005965 

526.0964355 

0*1327353 

-0*0002843 

535.0964355 

0*1352539 

0*0004914 

549* 0964355 

0*1379395 

-0*0003089 

556.0964355 

0*1411133 

0*0011219 

563.0564355 

0*1406250 

-0*0011093 

570.0964355 

0*1416016 

-0*0018787 

577.0964355 

o«i4379ee 

-0*0014214 

5 84 . C964355 

0*1469727 

0*0000095 

591.0964365 

o*i472i oe 

-0*0014893 

596.0964355 

0*1481934 

-0*0022557 

608. 0 964355 

0*1513672 

-0*0006246 

612.0964355 

0*1516113 

-0*0023236 

619.0964355 

0*1560059 

0*0003279 

626.0964355 

0*1574249 

0*0000041 

633.0964355 

0*1583862 

-0*0007776 

649.0964355 

0*1 604767 

-0*0004300 

647.0064355 

0*1656189 

0*0029692 

654.0964355 

0*1679230 

0*0035303 


DERIVATIVE OF THE ARGUMENT OF PE 
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Graph 1.- -Derivative of the argument of perigee as a function of the 
number of days past epoch. 
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Graph 2. — Computed argument of perigee as a function of the number of 
da> s past epoch. 
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Graph 3. — Derivative of the longitude of ascending node as a function of 
the number of days past epoch. 
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Graph 5. --Difference between the computed and observed values of the 
argument of perigee as a function of the number of days past epoch. 


COMPUTED OBSERVED 


iiiiiiiiiiiiiiiiimimiiiiiiiiiiliiiiliilillijl 


iiiiiimiiiiiiiimiiiiiiimiiHiiiiiiiiiiiiiiiiiiiisii 

iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiUiiriiiiii 

iiiiiiiiiimiiiimimiiiiiiiimiiiiiiiiiiKiriiiiiiii 
uumuuuiuuumiuuuimmuumiwiuium 

muuuuummmuuuumuuuumHiuuium 
iiiiiiiimiiiiiiiiiiiiiiiiiiiiiiimiiriiiiiiiiiiiiiiiii 

iiiiiumimiiuuiimiumiiii'miiumiiiiiiiim 
: i ii i 1 1 1 i i i imumiumillllHl’ jlllllllllllll i a 1 1 1 1 1 1 1 
luuuummmummuuumuumimmifiuui 
iiiiiiiiiiiiiiiiiiiiimiiimziiiiiiiiiiiiiiniimiiiii 

IIIIIIIIIIIIIIIIIIIIIIIIIIIFMIIIIIIIIIIIIIIIIIIIIIIIIIII 

luummmuuuuuumiumummuuuumiiil 

lllllllllllllllllllllll>Sfllllllllllllllllllllllllllllllllllil 
lllllllllllllllllllll'i 




■■■imiiiiiivmmmmmmmmammmmi 

[ i u m i iiiiiiiiFfliiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii II 


Iiiiiiiiiiiiiiiiimiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii 


muuuurmiuummiimumiuuummuuim 

miiiiunummiumiummimmumiuimim 

mmi’:uumummuuumuuuuumiuuium 

murjmummimmummmmiiimimimm 


mu 

mu 

mu 

mu 

■■■■I 

mil 

mil 

8 S 8 I 8 

mil 

mil 

mu 

mu 

iiiii 

mu 

mu 

mu 

mu 

mu 

mu 

urn 

urn 


mu 

mu 

iiiii 


luriiuumumiuumumuiM 
luuuuuuiuuuimuuummiil 
Piiuuiiiimiiiiiiiiiiiiiuimiiiiiii 


mi 


iiimui 
■mm 
iiiiiii 


100 200 

31 


NUMBER C 

fference bet) 

,veen the 

lg node as a 

function 


iiiii 
iiiii 
iiiii 
iiiii 
iiiii 
iiiii 
iiiii 
iiiii 
iiiii 
mu 
mu 
mu 
mu 
iiiii 
iiiii 
■iiiii 
uuiiuui 
uuiiuui 

UUIHIIII 

UUIIUUI 


Ml 


e nu 







PERIODIC PART OF (9 computed ~ 9 OBSERVED' (Degrees) 


1 .> 


. 0 ) 4 o 

.0125 

.0104 

.0083 

.0062 

.0042 

.0021 

0 

-.0021 

-.0042 

-.0062 

-.0083 

-.0104 

-.0125 



Graph 7. ---Periodic part of the differonce between the computed and ob- 
served values for the argument of perigee as a function of time. 
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Graph 8. --Least squares fit to the periodic part of the difference between 
the computed and observed values for the crgument of perigee as a function of time. 
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Graph 10. --Periodic part of the difference between the computed and ob- 
served values for the longitude of ascending node as a function of time. 
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